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A formulation is developed for general relativistic ideal magnetohydrodynamics in stationary 
axisymmetric spacetimes. We reduce basic equations to a single second-order partial differential 
equation, the so-called Grad-Shafranov (GS) equation. Our formulation is most general in the sense 
that it is applicable even when a stationary axisymmetric spacetime is noncircular, that is, even when 
it is impossible to foliate a spacetime with two orthogonal families of two-surfaces. The GS equation 
for noncircular spacetimes is crucial for the study of relativistic stars with a toroidal magnetic field 
or meridional flow, such as magnetars, since the existence of a toroidal field or meridional fiow 
violates the circularity of a spacetime. We also derive the wind equation in noncircular spacetimes, 
and discuss various limits of the GS equation. 

' PACS numbers: 04.20.Cv,04.40.Dg,52.30.Cv,95.30.Qd,95.30.Sf,97.10.Cv 

o : 
o 

^ , While most neutron stars have magnetic fields of ~ 10^^-10^^ G, studies of soft gamma-ray repeaters and anomalous 
' X-ray pulsars suggest that a significant fraction (> 10%) of neutron stars is born with larger magnetic fields ~ 10^''- 
The internal magnetic field of a new born neutron star may be even larger > 10^^ G if it is generated 
by the helical dynamo 0,0- A magnetic field of nearly maximum strength allowed bythe virial theorem ~ 10^^ G 
' may be achieved if the central engine of gamma-ray bursts are neutron stars 0, IMlS • In such super- magnetized 
1 neutron stars, so-called magnetars the magnetic fields have substantial effects on their internal stellar structure. 

■ Especially the deformation due to the magnetic stress becomes non-neglig ibie naiiiiiiiii. Since the deformation 

affects the precession, oscillations and the gravitational wave emission of neutron stars [ill IT^ [T3. Hsl IT^ . it is 
CNJ = important to investigate equilibrium configurations of magnetars. 

I General relativistic effects are sizable in the interior of a neutron star, so that any quantitative investigation of the 
• magnetars has to be based on general relativistic magnetohydrodynamics (MHD) p^Tl lla . [l9ll20| . Therefore, wc have 
<^ [ to solve the matter and electromagnetic field configurations in a curved spacetime, and have to take account of the 
O ■ electromagnetic energy-momentum as a source of the gravitational field. So far several works have been devoted to 
^ '■ equilibrium configurations of a magnetized star in a stationary axisymmetric spacetime ^3 1^ 113 ■ However these 
$H . works consider only poloidal magnetic fields for simplicity, since the existence of only a poloidal field is compatible 
with the circularity of the spacetime |23ll24| . In a circular spacetime, there exists a family of two-surfaces everywhere 
orthogonal to the plane defined by the two Killing vectors associated with stationarity 77'' — {d/dt)'^ and axisymmetry 
— [d/dipY |2^|2^|23|- Thus one may choose the coordinates {x^) — (t, x^,x'^,Lp) such that the metric components 
501; ff02, ff3i and 332 are identically zero. As a consequence, the problem is simplified dramatically. 



. However non-negligible toroidal magnetic fields are likely to exist in nature. Differential rotation generated during 
■ " ■ the gravitational collapse [28L or in the binary coalescence [s^l may wind up the frozen-in magnetic field to 
amplify the toroidal component 0, |3 • A toroidal magnetic field may be generated by the a-f2 dynamo during the 
first few seconds after the formation of a millisecond pulsar |2b tSj • In addition, convective motion may also exist in the 
interior of a neutron star [sj], which also violates the circularity of the spacetime [2^|3- Thus, we have to consider 
noncircular spacetimes. The degree of noncircularity of the spacetime in a neutron star with mass M^, and radius i?,, 
will be about ^ {M^/ R^)vmi and ~ 0.1-0.01(M,/i?*)7?.M, where v^i is the velocity of the meridional flow and TZm is 
the ratio of the magnetic energy to the gravitational energy [s^ • 

The problem to obtain an equilibrium configuration of a magnetar can be separated into two parts. The first 
part is the Einstein equations which determine the spacetime geometry under a given configuration of matter and 
electromagnetic fields. The second part is the matter and electromagnetic field equations in a given spacetime 
geometry. A (2 -I- 1) -I- 1 formalism to solve the Einstein equations under the presence of a spatial Killing vector was 
developed by Maeda et al. and by Sasaki [s^. This formalism is similar to the well-known 3-1-1 formalism [4^ . 

Later Gourgoulhon and Bonazzola '23'| developed a similar but different (2 -I- 1) -I- 1 formalism which is more suited 
for stationary axisymmetric spacetimes. So here we focus our attention on the second problem, i.e., to formulate the 
equations of motion of matter and electromagnetic fields in a curved spacetime. 

It is well known that the basic equations for a stationary axisymmetric ideal MHD system can be reduced to a 
single second-order, nonlinear partial differential equation, the so-called Grad-Shafranov (GS) equation, for a quantity 
called the flux function, ^. The GS equation was derived in the Newtonian case the Schwarzschild spacetime 
case (3^, and the Kerr spacetime case [s^ [s^. The flux function 5* is such that it is constant over each surface 
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generated by rotating the magnetic field lines (or equivalently the flow lines) about the axis of symmetry and the GS 
equation determines the transfield equilibrium. Any physical quantities can be calculated from the solution ^' of the 
GS equation. However, the GS equation in noncircular spacetimes has never been derived explicitly. 

In this paper we derive the GS equation explicitly in noncircular (i.e., the most general) stationary axisymmetric 
spacetimes. This is a first step toward the study of equilibrium configurations of magnetars. An attempt to solve the 
GS equation will be discussed in a subsequent paper. This paper is organized as follows. In Sec . HH we briefly review 
the conservation laws in stationary axisymmetric general relativistic ideal MHD systems '201 that are used to 
characterize the matter and electromagnetic field configurations. We neglect dissipative effects, which is a reasonable 
assumption because of the high conductivity and the low viscosity in neutron stars. In Sec. IIIII we derive the GS 
equation in an un-elucidated form. At this stage it is not clear if the GS equation is a second-order differential equation 
for the flux function ^. In Sec. lIVI we briefly review the (2 + 1) + 1 formalism by Gourgoulhon and Bonazzola to 
describe the geometry of noncircular stationary axisymmetric spacetimes in a transparent way. We do not, however, 
discuss the Einstein equations but assume the geometry to be given. In Sec. we explicitly demonstrate that all 
physical quantities except for the metric can be evaluated from the flux function ^E". In Sec. I VII we write down the 
GS equation in the covariant form projected onto the 2-surface T,tip spanned by t =const. and ip =const.. We also 
discuss various limits of the GS equation. Finally, we summarize our result in Sec. lVIll The energy-momentum tensor 
decomposed in the (2 + 1) + 1 form is given in Appendix^ and notation and symbols are summarized in Appendix IbI 

We use the units c — G — ks — ^- Greek indices (/i, v, a, P, • • • ) run from to 3, small Latin indices (z, j, k, • ■ ■) 
from 1 to 3, and capital Latin indices {A, B, C, ■ ■ ■) from 1 to 2, where — t and = The signature of the 
4-metric is (— , +, +, +). 

II. BASIC EQUATIONS AND CONSERVATION LAWS 
A. Basic equations for general relativistic magnetohydrodynamics 

The basic equations governing a general relativistic ideal MHD system are as follows [13, m. Baryons are conserved, 

(pu^),^^u^'p^,,+ pur^^, = 0, (2.1) 

where p is the rest mass density (i.e., the baryon mass times the baryon number density) and is the fluid 4- velocity 
with 

Uf,u^' = -1. (2.2) 
The electromagnetic field is governed by the Maxwell equations, 

= 0, (2.3) 

F^'".,, = 47rJ'', (2.4) 

where F^i, and are the field strength tensor and the electric current 4- vector, respectively. Equation ()2.3|) implies 
the existence of the vector potential A^,, 

The electric and magnetic fields in the fluid rest frame are defined as 

E^, - F^,u\ (2.6) 
= -ie^'-^^u.i^a^g, (2.7) 

where e^^ap is the Levi-Civita antisymmetric tensor with 69123 = V^d- Equations H2.6|l and H2.7|l are inverted to give 

F^, = u,,E, ~ u,E^ + t^,^0u"B^, (2.8) 
with EfjtU^ = Bfj,u^ — 0. In the ideal MHD, we assume the perfect conductivity, so that 

E^, - F^^u'^^Q. (2.9) 
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The equations of motion for the fluid are given by T^^".^ = 0, where T''^ is the total energy-momentum tensor of the 
fluid and electromagnetic fields, 

T^"^ = ip + pe+ p)u^u'' + pg^''' + - \g>"'F''^F^^ . (2.10) 

Then we obtain the conservation of the fluid energy, 

u^'{p + pt),^ + {p + pe + p)u^',^ = 0, (2.11) 

and the Euler equations, 

{p + pe+p)u^';,u'' + [g^"" + u^u'')p.,. - F^" J, = 0, (2.12) 

where e and p are the internal energy per unit mass and pressure, respectively. Eliminating u^;^ from Eqs. H2.1|l and 
(jnT|) gives 

u'^(p + pe),^ ^Aiu'^p,^, (2.13) 

where 

Ai = l + e+-, (2.14) 
P 

is the enthalpy per unit mass. Assuming local thermodynamic equilibrium, the first law of thermodynamics is given 
by 

de = -pd[^ +TdS, (2.15) 

where S and T are the entropy per unit mass and the temperature. Then Eqs. (|2.14(l and l|2.15|l imply 

dii^ — +TdS. (2.16) 
P 

Finally we supply the equation of state, 

P^pip,S). (2.17) 



B. Conservation laws in a stationary axisymmetric spacetime 

Here we recapitulate the conservation laws derived from the basic equations in the previous subsection in a stationary 
axisymmetric spacetime. There exists two Killing vectors associated with stationarity and axisymmetry, which we 
denote by rj'^ and respectively. The Lie derivatives of all physical quantities along the Killing vectors must vanish, 
e.g., C^u^^ — ^"^u^-u — u'^^^-i, = 0. We take r]f^ — {d/dty and = {d/d(pY so that x° — t and — (p are the time 
and azimuthal coordinates associated with the Killing vectors rj^ and respectively. Thus all physical quantities 
are independent of t and i£. 

Bekenstein and Oron [l9l l20j showed that a stationary axisymmetric system has several conserve d q uantities 
along each flow line. This is a general relativistic generalization of Ferraro's integrability condition |43,|43,|43- By 
exploiting the gauge freedom to make A^^^rf = A^^ = and ^ A^^a — 0, we can show that the magnetic 

potential 4" := A^£^^ = A3 as well as the electric potential $ := ^1/^??^ — Aq are constant along each flow line, i.e., 
u'^{S^'^ Ai^)^fj, = u'^{r]'^ A^)^fj, = Henceforth we label the flow line by which we refer to as the flux function as 
in the non-relativistic case [s^. The = const, surfaces are called the flux surfaces, which are generated by rotating 
the magnetic field lines (or the flow li nes) about the axis of symmetry. 

According to Bekenstein and Oron [l9l |20| , one can show that 



^03 


= 0, 


(2.18) 


FoA 




(2.19) 


F31 


= = C^f^pu^, 


(2.20) 


F23 


= *,2 = C^f^pv}, 


(2.21) 


F12 




(2.22) 
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where ft{^) and C{^) are conserved along each flow Hne and hence are functions of the flux function Vf. The above 
equations are effectively first integrals of the Maxwell equations H2.3|l . It may be useful to rewrite the above equations 
as 

= -Cp [{uo + nu3)u'' + + n^"] . (2.23) 

Note that ri(^') is the ^P-derivative of the electric potential, = —d^/d'^, and coincides with the angular velocity 

d(p/dt = /vP = il if there is no toroidal field F12 = and C 7^ 0, from Eq. (|2.22|) . In addition, one can show that 
i?(^), L{'^) and £'(^') are also conserved along each flow line 19], where 

-D = p{uo + nu3), (2.24) 

-E - (^f,+ ^^uo + C{uo + nu3)^, (2.25) 

L - (^^i+^^U3 + C{uo + nu3)^, (2.26) 

and B^^B''. Equations (|2.24|) . (|2.25|) and (|2.26|l are first integrals of the equations of motion T'^'';„ = along 

_Bp, 7]^ and respectively. In particular, Eq. (|2.25() is a generalized Bernouilli's equation. These conserved quantities 
are not mutually independent but there is a relation among them ^19j] . 

D^E-flL. (2.27) 

Together with Eqs. ifT^ ~ this imphes that 

B^ 

— + C(Ba + nBs) = 0. (2.28) 
P 

With Eq. we can rewrite £;(*) and /.(^f) in Eqs. TTI^i and as 

= /it^o-i^Cr!(uoS3-U3So), (2.29) 
47r 

1 



AiU3 + -^C(uoS3-M3Bo). (2.30) 

47r 



Finally, from Eqs. H2.13|l - (|2.15() one finds that the entropy per unit mass S is conserved along each fiow line, 
u^S,^ = 0, as a result of the perfect fluid form of the energy-momentum tensor. For a stationary axisymmctric 
spacetime, this implies that S is a function of 5 = S{^). 

In summary, for a give flux function ^, there exists flve conserved quantities, i?(^), i(^), f^(^), C(^) and 
S'(^'). Except for S{^), there are no perfectly relevant physical interpretations of these quantities. Nevertheless, 
by considering several limiting cases, we may associate them with terms that describe their qualitative nature. We 
may call E{'^) the total energy per unit mass, L{^) the total angular momentum per unit mass, ri(^) the angular 
velocity, and C(5') the magnetic field strength relative to the magnitude of meridional flow. Since these conserved 
quantities are essentially the first integrals of the equations of motion, specification of these functions characterizes 
the configuration of the fiuid flow and the electromagnetic field. As we will see in Sec. all physical quantities are 
completely determined once these conserved quantities are given as functions of the flux function provided that 
the spatial configuration of ^' is known. Therefore the problem reduces to solving an equation for the flux function 
that determines the spatial configuration of 5", that is, the GS equation. 



III. GRAD-SHAFRANOV EQUATION IN THE COMPONENT EXPRESSION 

The GS equation is given by the transfield component of the Euler equations (|2.12() . In this section, we consider 
the a;"^-derivative of the flux function in the Euler equations, and factorize the resulting equation to derive the GS 
equation. We express equations in terms of their explicit coordinate components, since it is the most straightforward 
way to incorporate the symmetry, e.g., (• ■ • ),o = (• • ■ ),3 = 0. Accordingly, the GS equation is given in the component 
expression. A covariant form of the GS equation based on (2 + 1) + 1 formalism, which may be useful for numerical 
calculations, will be given in Sec. I VII 

Using Eq. H2.16|l . the Euler equations H2.12(l can be written as 

Pfiu^-^u" + pfi^f^ + pUf^u^fj.^^ - F^^^J" - pTS^^ = 0, (3.1) 
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where we have used that u'^S'.^ = 0. First, let us consider the x^-component of the Euler equations (|3.1|) . The first 
term can be expressed as 



where FJ^^ is the Christoffcl symboL With Eq. (|2.2|) . the third term in Eq. H3.1|l can be transformed as 

— —pp^l{l + Uou" + U2U^ + UzU^) + puiu'^p^2 

= pU^{uxP.2 - U2P,l) - PP,1 - pU^{pUo)s - pU^{pU3)^l 

+ppu°{uo.i + F(,\wp) + ppu^{uz-i + F^s^Wp) 

= pu^ [{pUi).2 - {p-U2)-l)] - pp.^1 - pu"{pUQ)^i ~ pU^{pU3)^i 
-pp{u^Ui.i + U^Ui.2) + pi^iu^V^iU^ + U^F^'gMp), 

where we have used the fact u^^u^-^i, = 0. Thus putting Eqs. (|3.2|) and H3.3|) together, and using Eq. 

p{p.Ul-^U^ + ^,1 + WlU^/i,^) = pu^ [{pUl)-2 - (M"2);i] - pu° (pUo) - pU^{pU3)^i, 

[(M'"i);2 - (MU2);i] ^',1 - pW°(/iMo),l - pU^ {pUs) ^ 



Next, the fourth term in Eq. H3.1|l can be transformed as 



(3.2) 



(3.3) 
we find 

(3.4) 



- j3) 



AB-, 



Fi2iV^F'').i, 



where the second hue follows from Eqs. (|2.19|l and H2.2()(l . From Eqs. (|2.8|l and H2.9|l . we have i^^^ = 
U3B0)/ ^y—g. Therefore, together with Eq. (|2.22|l that gives i^i2, the last term in Eq. I|3.5|l is expressed as 

-^^F,2{V^F'^)^l = --^Cp(u=^-17u°)(7/oS3-U3So),l 

= --LpMO(?/oS3 - U3So)(C17)'*,i + -Lpu3(7ioB3 - UsBo)C'^,l, 

An in 



(3.5) 

-(U0-B3 - 



pu 



An 



pu 



■j-C{uoB3 - U3B0) 
An 



where primes denote differentiation with respect to Finally, the last term in Eq. (|3.1f) gives 

-pT5i = -pTS'-^^i. 

Combining Eqs. H3.4|l - (|3.7|l . we have 
f 1 



[ipui)^2 - {^^u2)s] - (J' - J°) 



An An 



'W,i 



(3.6) 



(3.7) 



-pu 



puo - -j-Cil {U0B3 - U3B0) 
An 



pu 



pu-i + -5-C (U0S3 - ^3-60) 
An 



0. 



(3.8) 



Recalling the expressions for E and L given by Eqs. I|2.29|l and l|2.3()|l . respectively, we see that the last two terms are 
just their derivatives. Therefore, we can factor out the a;^-derivative of the flux function in Eq. (|3.8|l to obtain 



[(/iui),2 - {pu2),i] - - nj°) 



-pu 



E'- — {U0B3 - U3B0) {cny 

An 



pu 



L'~ — (uoBs - U3B0) C 
An 



pTS' = 0. 



(3.9) 
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The same analysis applies to the a;^-component of Eq. I|3.1() , and one finds the above equation (|3.9() with the replace- 
ment of 1 by 'i' 2- Therefore, by assuming ^ ^ Q {A = 1, 2), the GS equation is given by 

- nj° + --^ iifiui) 2 - {piU2) i] - pu° \E' - A (Cn)'] + pu^ [L' - AC] + pTS' = , (3.10) 

where, for convenience, we have introduced an auxiliary quantity A defined by 

K= ^{uoB^-mBo). (3.11) 
47r 

At this stage, however, it is not clear if Eq. H3.10|l gives a second-order, nonlinear partial differential equation for the 
fiux function since the dependence on the fiux function is unknown. In Sec. [3 we explicitly demonstrate that all 
the physical quantities appearing in the above equation can be expressed in terms of 4' and its x'^-derivatives, and 
in Sec. I VII we derive the GS equation in the covariant form and we make it explicit that it is indeed a second-order, 
non-linear differential equation for 



IV. (2 + 1) + 1 DECOMPOSITION 

In this section, we briefiy review the (2 + 1) + 1 formalism of the Einstein equations for stationary axisymmetric 
spacetimes developed by Gourgoulhon and Bonazzola 23] , in order to describe our metric in a covariant fashion. Note 
that this formalism is different from the (2 -I- 1) -I- 1 formalism by Maeda, Sasaki, Nakamura, and Miyama j33, 34] and 
Sasaki (35|. which is suitable to the axisymmetric gravitational collapse. Here we adopt the formalism by Gourgoulhon 
and Bonazzola because it is more convenient for a spacetime which is not only axisymmetric but also stationary. 

Let be the unit timelike 4-vector orthogonal to the t — const, hypersurface Et and oriented in the direction of 
increasing t, 

= -iVi,p. (4.1) 
The lapse function N is determined by the requirement 

r^n^ = -1. (4.2) 

The 3-metric induced by g^n on Et is given by 

h^u = g^iv + n^_,ni,. (4.3) 

Similarly, let be the unit spacelike 4-vector orthogonal to the t = const, and Lp = const. 2-surface Yit^ and 
oriented in the direction of increasing (p^ 

= Mh^^ip^^ = Mip^f,, (4.4) 

where the vertical stroke ] denotes the covariant derivative associated with the 3-metric h^,y. The coefficient Af is 
determined by 



1. (4.5) 



The induced 2-metric on Y^np is given by 

H^^iy = h^^ - m^m^ = g^j_y + n^Ui, - m^m^. (4.6) 

The covariant derivative associated with the 2-metric H^i, is denoted by a double vertical stroke ]|. There is a relation 
between the determinants as 

= NVh = NMVh. (4.7) 

Any 4-vector can be decomposed into its projection onto T,t,p, the component parallel to and that to m^. The 
Killing vectors are decomposed as 

ri^' = Nn^" - TV'' = Nn^" - MN'^m^' - iVs'', (4.8) 

= Mm>'-MT.^', (4.9) 
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where the shift vector is (minus) the projection of ij'^ onto St, M^^ is (minus) the projection of f onto Et^, 
and N's'^ is the projection of N'^ onto 'Sfip- For our choice of coordinates, i.e., for =t and — 93, the component 
expressions for n'^, n^, and rn^ are 



(-7V, 0,0,0), 

1 m 

(-M7V^,0,0,M), 
Ms^ Ms^ 1 



0, 



M 



M ' M 



Note that N^f" = (O, TVe^, TVs^ O) and = TVs^ + N^^M^^. 

The exphcit component expressions of g^", g^i,, /i^", and /i^i, are given by 



9m goj 

9i0 9ij 



hAB hA3 



^33 ^33 




hij 

A T\,f_B 11 ,f_ A 



Af2 



V 



My 



M2 



M2 
1 

M2 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



where k — 1,2,3 and A,B — 1, 2. We can express the 4-metric g^^, in terms of N, iV^, Nj^"^, M, M^'^ and Hab 



- 2N^Adtdx^ - 2My:Adipdx^ + HABdx^dx^ + (M^ + AfsAMs^^) 



(4.18) 



where the functions N, N^, N^^, M, and iJ^is depend only on the coordinate {x^,x^). Since we only assume 

that physical quantities are independent oi x'^ = t and x'^ = (p, the metric in Eq. (|4.18l) has some freedom in the 
choice of coordinates. We will leave the coordinate freedom unspecified. In Scc. lVII the covariant GS equation will be 
given as an equation projected onto Eii^. 



V. PHYSICAL QUANTITIES FROM FLUX FUNCTION * 

Provided that the metric g^^ is given and the conserved quantities -E(vE'), i(^'), 5^(^), C(^') and S{'^) are given as 
functions of all the physical quantities can be evaluated once the (effectively 2-dimensional) configuration of the 
flux function * is known (see [i^ for the circular case). In this section, we explicitly demonstrate this fact for the 
most general case of noncircular spacetimes. 
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A. Fluid 4-velocity u'^ 

First let us consider the fluid 4-velocity . It is useful to prepare two vectors 77'' + fi^'' and + Qrj^^ constructed 
from two Killing vectors 77^ and and to make them orthogonal to each other (77'' + + Bry^) = by taking 

Then we can decompose the fluid 4-velocity in the coordinate bases as 

u^' = w„(77^ + r!e^) + u^{i^' + 677^) + US'", (5.2 

where 77'^ + fi^^ = (1, 0, 0, fi), + 877^ = (9, 0, 0, 1) and — (0, u^, u^, 0) in the component expressions, and hence 
UT,^nf^ = u^i^'-m^ = from Eqs. (|4.1U|I and l|4.12|l . and 

u° = u,, + Oii^, (5.3 
= + Qufj. (5.4 

The decomposition in Eq. (|5.2|) is not conforming to the spirit of the (2 + 1) + 1 formalism but makes it easy to obtain 
the coefficients and as shown below. 

From Eqs. H2.20|l and (|2.21|l . the term is given by 

where the antisymmetric tensor e'"' is defined by 

e'"' = e'"'"'^n„mi3. (5.6 
With Eqs. (|2.23|) - H2.26|l . (|4.8|l and (|4.9I) . the coefficients and are expressed as 

E-flL Ny. D Ny. 



where G,, and are defined by 



and Ny, and My are defined by 



Note that, using Eq. (|5.5|) . iVs and My are expressed as 

1 



Note also that 



Af 2 . 47r7i 



(5.7: 



_ _{J^-_QEl( J^\( \My 



G., = -(7/^ + ^!ep)(7/'' + ^7e^) = -(goo + 2^!go3 + ^^'533), (5.9 
Q = (Cm + 0^m)(C^ + =533 + 29503 + 6^500, (5.10 



Ny = Uy^{Ny^+nMy^), (5.11 
Ms = Uyf'iMy^ + eNy^). (5.12 



A^s = + r!MsA) , (5.13 

= II a (Ms A + ONyA) . (5.14 



GriG^p 

^ f^Geue^ + Tis^TisA - 2u^My + , (5.15) 
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is the square of the effective Alfven Mach number MAif (where the second equahty follows from Eg. (|5.25(l below). 
At the Alfven point AfAif = 1, the numerator L — QE should vanish to keep the velocity finite |4g. 

Thus, from Eqs. (|5.2|l . (|5.7|l . (|5.8|l and H5.5(l . given the metric 5^1, and the conserved functions L{^), ^{^), 

C(5') and the fluid 4- velocity can be obtained from the flux function 5" and its first derivatives _a if the 

density p and the enthalpy /i are additionally known. The expression for p will be given in the next subsection. The 
enthalpy p is then determined as a function of p and ^' as will be discussed also in the next subsection. 

Once the components of are known, the (2 + 1) + 1 decomposition of the fluid 4-velocity is easily performed. 
With the help of Eqs. and we have 



UnW' 



where 



(5.16) 



M [{n - N'P) + (1 - iV^e) uj] 



u^'' - (u„ + eu^)N^' 



{u^ + Q,u^)My_ 



(5.17) 
(5.18) 
(5.19) 



B. Density p and other thermodynamical quantities p, e, p and T 

The pressure p, the internal energy e, the enthalpy p and the temperature T are functions of the density p and the 
entropy S from Eqs. (|2.14() . (|2.15() and H2.17|l . Hence, given as a function of 5", the only remaining quantity to be 
known is the density p. 

The density p is determined by the normalization of the 4-velocity u^^u^ — ^G,fUrj^ + G^u^^ + u^aus"^ + 
2urjUs'' + + SujMs'' (Cm + ©^m) = -1' t^^at is. 



G„m2 ^ G^GYC^p^ y"G^p) ^ N^PC^p^ ' G 



B , iVs^ Ms' 



Gf 



(5.20) 



This equation is what is called the wind equation (see |45| for the circular case). Note that p contains the first-order 
derivatives through this equation. 

C. Magnetic field B^' 

The magnetic field is also calculated from the flux function. With Eqs. H2.23|l . H5.16() . (|4.8|l and (|4.9|l . the (2 + 1) + 1 
decomposition of the magnetic field is given by 



where 



B„ = CpN ^[G^u.^ + N^'jiu.r^ + Qu^)-! 
Bm = CpM (g^u^ + N^^ {{n- N'^)u,^ + {1- N'^e)u^} + N'P -n 



The magnetic strength is given by 



^2 ^ ^Mg^ ^ (72^2 



(g„u^ + ^s) ' - G„ = GV' [g„ (Geu^2 ^ u^^u^A _ 2u^M^^ + iVs' 



(5.21) 

(5.22) 
(5.23) 
(5.24) 

(5.25) 



The (2 + 1) + 1 decomposition of the energy-momentum tensor is given in Appendix IXI 
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D. Electric current 



Let us consider the following components of the electric current, 



I f 



AnNMVH 



NMVHF' 



1 



AttNMVh ^ 



A AnNM 
1 

,A ~ inNM 



{NMF' 



0A\ 



\\A' 



\\A- 



(5.26) 
(5.27) 



The field strength tensor components appearing in the above equations are also expressed in terms of the flux function 

as 

p.A ^ ^gOOgAB _ gOBgAO^ ^^^^ ^ ^^^Oi^^^S ^ ^03^^^^ ^ ^g01gA2 _ ^02^^!^ 



1 



H 



.4B , MiMi 



F 



3A 



1 



£B 



^3B A3 „33 „AB\ 



M2 



Fl2 

71' 



^g^OgAB „ g3BgA0-^ ^^^^ ^ (^SiS^^S „ g33gAB-^ ^ ^g31gA2 _ ^32^^!^ 



1 



where F12 is expressed as 



JP 



AB 



AP 

2\ 



Mi- 



iV2 



KP 



H 



AB 



7V2M2 
F12 

7h' 



= cNMp{u^ - f7u°) = CNMpu^ (1 - rje) . 



(5.28) 



(5.29) 



(5.30) 



In the above, the first equality follows from Eq. (|2.22(l . and the second from Eqs. H5.3|l and 1)5. 4|l . Thus, and 
are expressed in terms of ^ and its first and second derivatives. 



E. Auxiliary quantity A 

We also need to evaluate the auxiliary quantity A defined by Eq. (|3.11|) . that is. 



A= ^(uoBs-usBo). (5.31) 



From the expression of given by Eq. (|2.23|) . we have 

A = -T-C'P [uo (303 + ^^533) - -"3 (500 + ^903)] ■ (5.32) 
Using the component expressions of given by Eq. H5.2|l , this is rewritten as 

A = ~^'^P (^03^ ~ 500533) (1 - f^©) + ws^?7p (503 + f^533) - Ms^^M (500 + ^903)] 

= ^-^^P (.903^ - 500333) (1 - f^0) + A^s (500 + f^503) , (5.33) 

where the second line follows from Eqs. (|4.8|) . (|4.9II . H5.1|l and 15.12|l . The above form is sufhcient for A to be obtained 
from the flux function, but it can be further simplified if we use Eqs. H5.1|l and H5.10|l . From these equations, we find 

(1 - nQ) (500533 - .go3^) = Q (500 + f^503) • 

Therefore, we obtain 

A - ^ (tioSa - U3S0) - ^Cp {g^u^ ~ Afe) (500 + f^503) ■ (5.34) 



11 



VI. GRAD-SHAFRANOV EQUATION IN THE COVARIANT FORM 

Now we are ready to show that the GS equation H3.10|l is indeed a second-order differential equation for the flux 
function ^. At the same time, following the spirit of the (2 + 1) + 1 formalism, we express the GS equation in the 
covariant form with respect to the geometry of Sji^. 

The covariant expression for the GS equation is readily obtained as 

- nj" + ^^^e^''{fxu^A)\\B - p{u^ + [E' - AiCny] + p{u^ + nu,-,) [L' - AC] + pTS' = 0, (6.1) 

where we have replaced yj—g-, vP and in the original GS equation H3.10|l by their (2 + 1) + 1 type expressions (|4.7ll . 
(|5.3|) and H5.4|l . respectively, and, as before, a double vertical stroke |j denotes the covariant differentiation with respect 
to the 2-metric Hab- 

In the previous section, we have seen that u^, u^, us"^ and 9 (section [3 A), p, ^ and T (section [3 B), and 
(section 13 D), and A (section 0E) are all expressed in terms of ^ and its derivatives, given the conserved functions 
i?(^'), i('I'), f^('I'), C'(^) and S{'^), and the metric g^^^. In particular, we have seen that J° and contain the 
second-order derivatives of while p (hence p) as well as u^a contain the first-order derivatives of . Thus, the 
GS equation (|6.1|) is a second-order, non-linear differential equation for where the first three terms contain the 
second-order derivatives. 



A. No toroidal field limit 



From Eqs. H2.20|l - H2.22|l . we find that the toroidal field and the meridional flow vanish if |C| — > oo. Here note 
that - riu° = Mj(l - VLQ) cx C^^ ^ in Eq. (|2.22|l from Eqs. (|5.3|) . (|5.4|l and (|5.8|l (and hence coincides with 
the angular velocity d(f/dt — /vP = 0). In the absence of the toroidal field and the meridional flow, a spacetime is 
circular. The circular limit is expressed as 

iVs^ -> 0, Ms-^ -> 0. (6.2) 
Therefore, in the |C| cx) limit, the GS equation H6.1|l reduces to 

J3 - VtJ° - pu,, [E' - VLL' - CMl'] + pTS' = 0, (6.3) 
where the density is determined by {E — ilL)^ /Gjjp^ = G^u,,^ = 1 from Eq. H5.20|l . and 

Here we regard \C\ ^ oo as the limit of a sequence of models with \C\ < oo. The last term [l — (47r/i)/(G,,C'^/3)] ^ 
in Eq. (|6.4() can be neglected if the density p is finite. However, in the case when there is a surface with p — like 

a star and the flux function is not constant on that surface, the last term [l — (47r/i)/(Gr,G'^p)] diverges near 
the surface. Unless one can fine-tune the term L — QE so that its zero point cancels this divergence, which seems 
unlikely to be possible, we should demand the rigid rotation il' = in Eq. 1)6.3(1 . This is consistent with Bonazzola et 
ai. caiiii. Note that if the fiux function ^! is constant on the p = surface, we may find C that satisfies C p — + oo 
on the surface. 



B. No poloidal field limit 

The poloidal field vanishes if we let ^ 5'^ and take a limit 5 ^ 0, as we can see from Eqs. (|2.20|) and (|2.21ll . In 
this process we relabel the flow lines by 4' and replace the conserved quantities as ^ E{^), L L{'^), fl ^{^), 
C G(5') and S In the limit (5^0, the meridional flow vanishes us^ from Eq. 1(5.5(1 . Then we can 

show that the spacetime is circular as expressed in Eq. 1(6.2(1 [23| . Therefore, in the 5 — > limit, the GS equation 1(6. 1() 
reduces to 

{ur, + eu^) [e' - A {cn)'] - {u^ + nu,j) [V - ag'] - ts' = o, (6.5) 

where primes now denote differentiation with respect to VP. This is an algebraic equation. Here we regard \1/ — > as 
the limit of a sequence of models with ^' 7^ 0. If 5* is exactly zero, the transfield components of the Euler equations 
((2.12(1 are satisfied regardless of the GS equation (see Sec. IIIip . Therefore there may exist 'isolated' solutions which 
cannot be obtained by the limit discussed here. 
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C. No magnetic field limit 

There are two limits for configurations with no magnetic field. The first way to obtain such configurations is to let 
^ ^ (^i^' and C C/S2 and take the hmit ^ and 62 0, as we can see from Eqs. (|2.20f) - (|2.22|l . Here note 
that - fiu" = w^(l - 179) cx ^ in Eq. (|2.22|l from Eqs. (|5.3|) . (|5.4|l and (|5.8|l (and hence coincides with 
the angular velocity d(p/dt = /vP — Q). In this process we relabel the flow lines by ^ and replace the conserved 
quantities as £; -> E{^), L ^ L(^'), Q ^ r2(£), C £(£) and S S'(*). In the limit -> and S2 0, the 
meridional flow vanishes u^^ from Eqs. 12.2011 and 12.21|) . so that the spacetime becomes circular as in Eq. (|6.2|l . 
Therefore, the GS equation (|6.1|1 reduces to 

[{E - ilL)' - n' (CA - L)] ~ TS' = 0, (6.6) 

where primes denote differentiation with respect to ^, CA is given by Eq. (|6.4|l . and the density is determined from 
Eq. i^U^ as 



Gr,Ur,^ = 1 . (6.7) 



n ,,2 " 

Let us see the relation between this limit (5* — > and \C\ 00) and the case of a rotating star From 
Eq. I|6.7|l we have the Bernouilli's equation for a rotating fluid as 

In/i - Inu,, - ln(S - fiL) = 0. (6.8) 

From the definition of L, Eq. H2.30|l . we have CA — L = —^us = —u^Uri{E — ilL). Then for an isentropic star S' = 0, 
the GS equation in Eq. H6.6|l is written as 

U3Ur, = -^[~HE-nL)] , (6.9) 

where we have used the fact that — ln(_E — flL) can be regarded as a function of il, since E, L and fl arc functions 
of ^' only. Equation H6.9|l is the well-known integrability condition for a rotating fluid [21I . 

Note that if we regard |C| — > 00 as the limit of a sequence of models with |C| < 00, we should demand the rigid 
rotation ft' ~ 0, as discussed in Sec. lVI Al 

The second way to obtain configurations with no magnetic field is to let 4' S'^ and C — > SC and take the limit 
S ^ Q. In this process we relabel the flow lines by ^ and replace the conserved quantities as — > E{'^), L L{^), 
n ri(*), C (?(#) and S — > S{'^). As we can see from Eqs. (|2.20|) and (|2.21|l . there exists a meridional flow in 
this case. The toroidal field vanishes since in Eq. H2.22|l — ilu" = ^^(1 — flQ) and (L - QE)/G^^j, + M-^/G^ 

in the limit i5 ^ 0, from Eqs. (|5.3|) . H5.4|l and H5.8|l . The GS equation in Eq. (|6.1|l reduces to 

- pu,j [E' - flL') + pu^ [L' - BE') + ^-^£-43 (^usa)||b + pTS' = 0, (6.10) 

where primes denote differentiation with respect to 4*. One can introduce a yet new flux function ^ by d"^ — /G to 
absorb the function C into the definition of the new flux function ^ . The resulting equation may be directly obtained 
from the Euler equations for a perfect fluid. 



D. Newtonian limit 

In the Newtonian limit, all physical quantities are expanded in power series of the typical fluid velocity '4^ . The 
metric reduces to 

Of^^dx^'dx" = -(1 + 2<j))dt'^ + (1 - 2<j>) [dZ^ + dR^ + i?^ V) , 

where the Newtonian potential 4> is of order O(w^), and we adopt the cylindrical coordinate {t, Z, tp) 
the 3-dimensional velocity by 

J dx^ 
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where 

,,2 



= + y, (6.13) 



and v"^ = v^Vi. We regard the internal energy e and the pressure p to be 0{v^). To make the energy density of the 
electromagnetic field O(w^), we demand 

B'-O(w), B°r^O{v^), ^r^O{v), n^O{v), (6.14) 



from Eqs. and (|TTH|) - ifT^ . 
From Eqs. - we find 



= Cpi;^, (6.15) 
B"^ = Cp (v'^ - , (6.16) 

where B'^ := RB^ and v'^ := Rv^. From Eqs. (|On|) . ifO^ and lIT^ . we also have 

L = Rv"^ + CA = R(^v'f' - - 0(v), (6.17) 

2 

D-l = e+^ + ^ + (A-i?f7w'^-0(i;2). (6.18) 
P 2 

These results are to be compared with the Newtonian results (note the correspondences between our notation and 
that of [13 as f7 ^ G, C ^ An/F, L ^ -H/F, and D - 1 ^ J). 

Let us obtain the Newtonian GS equation, which is of order 0{v). First consider the and —flj'^ terms in 
Eq. H6.1|l . In the Newtonian order, these terms may be evaluated on the flat background with N — 1 and M = R^. 
Then — is found to be 0{v^), and the term is given by 

J3 = (—H'^^^b] = ^A**, (6.19) 



where 



did 92 

'''-^dRROR^OZ-^- (^-^O) 



Next consider the e (/iitsA)||s term in Eq. H6.1|) . By using Eqs. H5.19|l and IjS.Sfl . we have 



C^p C \Cp ' 



(6.21) 



Finally consider the terms proportional to p in Eq. H6.1|l . To the lowest order, from Eqs. I|5.1|l . H5.7|l . ()5.8|l and H5.34|l . 

we have 

e = R'^n^O{v), (6.22) 

= 0(1), (6.23) 
{L - R^n) / 47r \ / _ 47r 



= -^^^^j^l-^j -Oiv), (6.24) 

= (D - 1 + r^X)' ~ 0(u), (6.25) 

L' - 0(1), (6.26) 

A = -^^0(.). (6.27) 

Then we can show 

pu^ [E' - ACn' - nL'] = p[{D- ly + Rv^n'] . (6.28) 
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From Eqs. (|F?77|) and l|HT7jl . we may express p as 

47rA 



4tt L-Rv'p 

c 



Using this expression, we can show 



pu^ [V - AC] 



1 /47ri 



Rv^- 



47r\ C 



CR'^ \ C "'^ C } A-K 
1 ( 47rL ^47r 



/47rL 



47r 

-'If 



where the second equahty follows from Eq. Ht).17|l . It is also easy to show that pQu^ [E' — K{CVl)'] ^ 0{v'^). 
Therefore from Eqs. H6.19|l . (|6.21(l . H6.28|l and H6.30(l . the GS equation in the Newtonian limit is given by 



47r \ 



47r 

A*^' V 

C 



1 

C'p 



-inpR^ [{D - ly + Rv'^n'] 



UttL 



AnR^pTS'. 



(6.29) 



(6.30) 



(6.31) 



This is equivalent to the Newtonian GS equation in [36 



VII. SUMMARY 



We have derived the GS equation H6.1|l in noncircular (the most general) stationary axisymmetric spacetimes. The 
GS equation has been given in the covariant form projected onto the t =const. and if =const. 2-surface Et,^. We have 
also derived the wind equation (|5.2U|) in noncircular spacetimes. We have discussed various limits of the GS equation 
(no toroidal field limit, no poloidal field limit, no magnetic field limit and Newtonian limit). 

To obtain equilibrium configurations of magnetars, we have to solve the GS equation (|6.1|l . As first glance, it 
looks formidable to solve it. One possibility is to take a perturbative approach to solve the GS equation. Unless the 
magnetic field is as strong as the maximum magnetic field allowed by the virial theorem ~ 10^^ G |iri lll|. we may 
assume weak magnetic fields compared with gravity. Then the magnetic field may be treated as a small perturbation 
on an alr eady -known non-magnetized configuration. This approach is similar to that developed for slowly rotating 
stars |49l l5Clj |. in which the perturbation parameter is the angular velocity. Work in this direction is in progress. 
The preliminary study indicates that the degree of noncircularity of the spacetime in a neutron star with mass 
and radius R^, is about (iVs^A^s^)^/^ {M.^/R^)v^i and (Me^A/s^)^/^ - 0.1-0.01(A/*/i^*)7^A^ , where w^f is the 
velocity of the meridional flow, TZm is the ratio of the magnetic energy to the gravitational energy, and the length 
scale is normalized by the mass [32j . 
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APPENDIX A: (2 + 1) + 1 DECOMPOSITION OF THE ENERGY-MOMENTUM TENSOR 

The (2 + 1) + 1 decomposition of the energy-momentum tensor is 

T^"" = en'^n'' + j (n^m'' -I- m'^n'') + (n^ff^" + if^'^n'') 

-hsm'^m'' + {m^'HA'' + HA^m") + s^'^Ha^Hb" , (Al) 

where and are the unit timelike and spacelike normals to the 2-surface 'Stip, respectively, and A, B = 1,2. For the 
electromagnetic field in an ideal MHD system, we have from Eqs. and (ESJ, F^^F^'a = (u^'u" + g'^") B'^ - B^'B" 



and F'^'Ff,^ = 25^, and hence 



47r V 4" """J 47r 

Then, using Eqs (|5.1()|) and 15.21|l the components of the energy-momentum tensor in Eq. (|A1|) are obtained as 

1 



T^^n^^n^ ^ {p + pe+p) (ii„) -p + 



An 



(Un) 



2 1 1 n2 /n ^2 



1 



j = -T'^''nf_,m^ = {p + pe+p) UnU,n + ^ - B^B^n) , 

/ - -H^'^T^^'^n, = {p + pe+p) u^u^^ + i- {B^Ur^u^^ - B,,B^^) 



Tf^^m^m^ = {p + pe+p) (umf +P + 



47r 
1 



= H^^,T^"'m„ = {p + pe+p) UmU^^ + — [B^u.nU^^ - B^B^' 

in 



H^^H^'^T^^- = {p + pe + p) u^^u^'' + pH-' 



1 

An 



where u„, u™, us^, -Bn, -Bm, -Be'^ and B^ are given by Eqs. (|5.17|l - H5.19() and (|5.22|) - H5.25|l . respectively. 
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(A2) 

(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 



APPENDIX B: SYMBOLS 

Here, we summarize definitions of some of the symbols we use, which may not be commonly used, with the equation 
numbers where they are defined or introduced. 

• Quantities conserved along each flow line: 

• Flux function, * := A^^^ = A3. 

$ Electric potential, $ := Af^ij'^ = Aq. 

C 'Magnetic field strength' relative to the magnitude of meridional flow; Eqs. H2.20|l - ()2.22f) . 

D 'Fluid energy' per unit mass; Eqs. (|2.24|) and (|2.27|l . 

E 'Energy' per unit mass; Eqs. (|2.25() and H2.29|l . 

L 'Angular momentum' per unit mass; Eqs. H2.26|l and H2.30I) . 

il 'Angular velocity' of the magnetic field hne, fl — —d^/d^; Eq. H2.19I) . 

S Entropy per unit mass; Eqs. (|2.15|) . 

• Quantities associated with the metric: 



Killing vector associated with stationarity, 77^ = (d/dt)'^; Eq. (|4.8|l . 
Killing vector associated with axisymmetry, — {d/dtfY\ Eq. H4.9|l . 
Unit timelike 4- vector orthogonal to t =const. hypersurface Ej; Eqs. H4.1|l . 



h^v 3-metric on Sf; Eqs. (|4.3|) . 



Gri 

e 



Unit spacelike 4-vector orthogonal to t =const. and ip =const. hypersurface St,^; Eqs. H4.4|l . 
2-metric on Sf^; Eq. H4.6|l . 
Norm of 77^ + f7^^; Eq. lf5l^ . 

Quantity such that + 677^ is orthogonal to 77^ + il^^; Eq. H5.1|l . 
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Norm of + 9r)^: Eg. ( jFlUll . 

N Lapse function, N = —rj^^n^; Eqs. H4.1|l . 
N^' Shift vector, iV^ = TVn^ - r;^; Eq. 
iVs^ Projection of A^^ onto Ef^; Eq. 
M = m^^^; Eqs. ES- 
A/e^ = Afm'' - Eq. 

• Quantities associated with the fluid: 
u'" Fluid 4- velocity; Eqs. 

p Rest mass density; Eqs. (|2.1UII . 

/i Enthalpy per imit mass; Eq. (|2.14(l . 

e Internal energy per unit mass; Eq. H2.10|l . 

p Pressure; Eq. (|2.1U|I . 

T Temperature; Eq. lfTT5|l . 

= (uO - eu^)/{l - ne); Eq. or Eqs. |Ol and 

= (u^ - f7u°)/(l - 176); Eq. or Eqs. ^ and (1^ . 

Un = -n^u^; Eqs. (|5.16(l and H5.17|l . 

= nif^u^^; Eqs. H5.16|l and (|5.18|) . 
us^ Projection of u'^ onto St^^; Eqs. (|5.16|l and H5.19|l . 
us^ (A = 1, 2) components of u^; Eq. H5.2|l . 

iVs A component of u^^ defined by Eq. H5.11|l . 
Ms A component of iij:^ defined by Eq. H5.12|l . 

• Quantities associated with the electromagnetic field: 
Ef^ Electric field in the fiuid rest frame; Eq. H2.6|l . 

Magnetic field in the fiuid rest frame; Eq. H2.7|) . 
Bn = -n^B^' = n^u,e^''"'5F„^/2; Eq. 
Brn = m^B^' = Ut,m,e^''''^PF^p/2■ Eq. (ESI). 
Bt.^" Projection of onto E*^, Eq. (EHU- 

Electromagnetic current 4- vector; Eq. H2.4(l . 

• Others: 

MAif Alfven Mach number; Eq. IjS.lSfl . 

A An auxiliary quantity defined by Eq. H3.11() . 
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